Search for CP Violation in the decays — Kstt'^ and — KsK^ 
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A high-statistics sample of photo-produced charm from the FOCUS (E831) experiment at Fermi- 
lab has been used to search for direct CP violation in the decays D"*" Kstt'^ and D"*" —> KsK'^ . 
We have measured the following asymmetry parameters relative to D"*" K~-k^-k'^: Acp{Ksi^^) = 
(-1.6±1.5±0.9)%, Acp{KsK+) = (-h6.9±6.0±1.5)% and Acp(KsK+) = (+7.1±6.1±1.2)% rel- 
ative to -D^ Kstt'^- The first errors quoted are statistical and the second are systematic. 
We have also measured the relative branching ratios and found: r(D+ — > K^n^) / T{D'^ — > 
K--K+TV+) = (30.60 ± 0.46 ± 0.32)%, r{D+ K°K+) / V{D+ K--R+-K+) = (6.04 ± 0.35 ± 
0.30)% and V(D+ ~* K° K+) / V{D+ ^ K^-k+) = (19.96 ± 1.19 ± 0.96)%. 

PACS numbers: 11.30.Er 13.20.Fc 14.40.Lb 



CP violation occurs when the decay rate of a parti- 
cle differs from that of its CP conjugate §. In the 
Kobayashi-Maskawa ansatz this arises due to the non- 
vanishing phase in the Cabibbo-Kobayashi-Maskawa ma- 
trix when the decay amplitude has contributions from 
at least two quark diagrams with differing weak phases. 
In addition final state interactions (FSI) must provide a 
strong phase shift. In the Standard Model direct CP vi- 
olation in the charm meson system is predicted to occur 
at the level of 10^"^ or below The mechanism usu- 
ally considered is the interference of the tree and penguin 
amplitudes in singly-Cabibbo suppressed (SCS) decays. 



In the decay Z3+ Ksn^ the Cabibbo favored (CF) 
and doubly-Cabibbo suppressed (DCS) amplitudes con- 
tribute coherently with, perhaps, a different weak phase.* 
In addition the isospin content of the DCS amplitude 
differs from that of the CF case so we can expect a 
non-trivial strong phase shift. Several authors have com- 
mented on the effect of K'^ mixing on the CP asymmetry 
for this decay mode and the possibility of using it to 



The charge conjugate state is impUed unless stated otherwise. 
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search for new physics ||, ^ . 

Differences in the non-leptonic decay ampHtudes of 
charmed mesons are almost certainly due to FSI. These 
effects tend to be amplified in the charmed system mak- 
ing it an ideal laboratory for their study jsj . The isospin 
amplitudes and phase shifts in D KK, D — *■ Ktt and 
D — > TTTT decays can be extracted from measurements of 
the branching fractions For example the magnitude 
of the 1=3/2 amplitude can be obtained directly from the 
D+ K°TT+ partial width §. 

Previous studies of D'^ i^s7r+ and KsK^ 
have concentrated on measuring relative branching ratios 
^. This paper reports the first measurement of the 
CP asymmetry for these decays. 

The data were collected during the 1996-1997 fixed 
target run at Fermilab. Bremsstrahlung of electrons and 
positrons with an endpoint energy of approximately 300 
GeV produces a photon beam. These beam photons in- 
teract in a segmented beryllium-oxide target and produce 
charmed particles. The average photon energy for events 
which satisfy our trigger is ~ 180 GeV. FOCUS uses 
an upgraded version of the E687 spectrometer which is 
described in detail elsewhere Charged decay prod- 
ucts are momentum analyzed by two oppositely polarized 
dipole magnets. Tracking is performed by a system of sil- 
icon vertex detectors in the target region and by multi- 
wire proportional chambers downstream of the interac- 
tion. Particle identification is performed by three thresh- 
old Cerenkov counters, two electromagnetic calorimeters, 
an hadronic calorimeter, and by a system of muon detec- 
tors. 

The £>+ K^n^TT^ decay is reconstructed using a 
candidate driven vertexing algorithm. A decay vertex 
is formed using the reconstructed tracks after which the 
momentum vector of the parent D meson is intersected 
with other tracks in the event to form a production ver- 
tex. The confidence level of the secondary vertex is re- 
quired to be greater than 1%. The likelihood for each 
charged particle to be an electron, pion, kaon or proton 
based on the light yield from each threshold Cerenkov 
counter is computed [Q. We demand that the Kaon 
hypothesis Wk, (i.e. — 2 ln(kaon_likelihood)), be favored 
over the pion hypothesis W-^ by AW = VF^ — Wk > 
1. We also make a pion consistency cut by finding 
the alternative minimum hypothesis Wmin and requiring 
Wmin — Wtt > —2 for both pions. We eliminate con- 
tamination due to Z?*+ — > K~TT~^)Tr~^ by asking 
that neither Ktt invariant mass combination lies within 
25 MeV/c^ of the nominal mass. 

The techniques used for Ks reconstruction are de- 
scribed elsewhere iQ. Because 90% of Ks decays occur 
after the Ks has passed through the silicon strip detec- 
tor we are unable to employ the same vertexing algorithm 
used to reconstruct the Z)+ K^tt^tt'^ decay. Instead 
we use the momentum information from the Ks decay 
and the silicon track of the charged daughter to form a 
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FIG. 1: Invariant mass plots for D'^ Kstt^ and 



candidate D vector. This vector is intersected with candi- 
date production vertices which are formed from two other 
silicon tracks. As a final check we force the D vector to 
originate at our production vertex and calculate the con- 
fidence level that it verticizes with the charged daughter. 
This confidence level must be greater than 2%. We re- 
quire that the momentum of the charged daughter be 
greater than 10 GeV/c, that the confidence level for it to 
be a muon be less than 1%, and that it traverse the entire 
length of the spectrometer. For the decay Kstt^ 
we demand that Wmin — Wtt > — 6 and W^r — Wk < 0, for 
— > KsK^ we ask that the kaon hypothesis be favored 
over both the proton and pion hypotheses by requiring 
Wp - Wk > and - Wk > 3. We remove electron 
contamination by ensuring that the charged D daughter 
links to only one silicon microstrip track. Electron pairs 
usually have a very small opening angle in the silicon and 
chamber tracks tend to link to both tracks. Checks for 
electron contamination of the Ks sample using the elec- 
tromagnetic calorimeters showed no significant effect. We 
use only Ks candidates which have a normalized mass^ 
within three standard deviations of the nominal value. 
Additionally, to reduce backgrounds in the KsK^ 
mode, we do not use the category of Ks decays which oc- 
cur downstream of the silicon where both Ks daughters 
lie outside the acceptance of the downstream magnet. We 
make the same cut on the Z?^ —t Kstt^ normalization 



^ The normalized mass is the difference between the measured and 
the nominal mass divided by the error on the measured mass. 
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FIG. 2: Invariant mass plot for KsK^ and 

D' KsK~ . The shaded area is the smoothed background 

shape from D+ ^ K*+ K° and _D+ ^ K'°K+. 
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FIG. 3: Invariant mass plot for K^tv^tv^ and 

D~ K'^TT~Tv~ . The shaded area is the third-degree 
polynomial background in the fit region. 



signal. When the Ks decays in the silicon detector we 
demand that all three tracks be inconsistent with origi- 
nating at the same vertex. This eliminates backgrounds 
from decays such as D'^ — > 7r~7r+7r+. 

For all modes we require that the production vertex 
have a confidence level greater than 1%, that the maxi- 
mum confidence level for a candidate-!? daughter track 
to form a vertex with tracks from the primary vertex 
be less than 20%, that the significance of separation of 
the production and decay vertices be greater than 7.5 
and that both vertices lie upstream of the first trigger 
counter. The momentum of the D must be greater than 
40 GeV/c. In Figures 0, |^ and || we show the invari- 
ant mass distributions for the decays A'^tt^, KsK^ and 
K~-K^'K'^ respectively. 

We construct the CP asymmetry, Acp as the difference 
in the yields, (corrected for efficiency and acceptance), 
of the decay in question divided by the sum. We must 
also account for differences in production between the 
£)+ and . To do this we ratio the corrected yields to 
those of a Cabibbo favored decay which is assumed to be 
CP conserving. We measure: 



A 



CP 



where for example. 



Tj{D+)-v{D-) 
7^{D+) + 7^{D-) 



N{D+ Kstt"^) 

N{D+ K-TT+TT+ 



is the ratio of the corrected yields for each decay which 
is equivalent to the relative branching ratio. 



To account for non-Gaussian tails in the — > 
K^TT^TT^ signals we find it necessary to fit these distribu- 
tions using two Gaussians and a third-degree polynomial. 

The KsTT'^ distribution is fit using a Gaussian and a 
linear polynomial. The non-linear background shape be- 
low 1.75 GeV/c^ in the Kstt~^ plot is primarily due to 
£)+ ^ Ksl'^vi and is not included in the fit. 

We fit the KsK^ signal using a combination of a 
Gaussian, linear polynomial, and a background shape 
derived from Monte Carlo. This shape is a smoothed 
fit to D+ K*+K° and D+ K*°K+ decays which, 
due to a missing 7r°, are responsible for the background 
shape below the D+ peak. Because of the difficulty in 
fitting the region between the D+ and Df peak in the 
KsK+ distribution we only fit up to 1.935 GeV/c^. To 
minimize systematic errors we change the Ks selection 
cuts on the D'^ — > Kstt^ normalization signal to match 
those used for the D'^ KsK^ mode. The yield for the 
decay Ks7t~^ changes to 4487±96 events and for 

KsTT^ becomes 4770±96 events. We can now cal- 
culate the relative branching ratios and CP asymmetries. 
The results are shown in Tables || and 

We studied systematic effects due to uncertainties in 
our Monte Carlo production model, reconstruction al- 
gorithm, and variations in our selection cuts. For the 
— + KsTT^ measurements we split the sample into 
eight statistically independent subsamples based on 
momentum, loose and tight normalized Ks mass cuts, 
and the time period in which the data were collected. 
The momentum dependence of the result arises mainly 
due to uncertainties in the parameters used to generate 
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TABLE I: Relative branching ratio results. The first error 
is statistical and the second is systematic. We account for 
the decay chain A'" — » Ks ■k'^-r~ by multiplying our 
Ks numbers by a factor of 2.91 assuming that T{D'^ — > 
_f5'°7r+) = 2 X r(_D+ KsTv^) ; we then quote these results 
in terms of A"". 



Measurement 



T(D 



Result PPG Average 



r(D+^K-.+.+ ) (30.60 ±0.46 ±0.32)% 

no'l^S/Z'+ll) (6.04 ±0.35 ±0.30)% 

r'olZf 0^+) (19.96 ± 1.19 ± 0.96)% 

"This is the measurement of reference ^ with statistical and sys- 
tematic errors added in quadrature. 



(32.0 ± 4.0)% 
(7.7 ±2.2)%" 
(26.3 ±3.5)% 



TABLE II: CP asymmetry measurements. The first error is 
statistical and the second is systematic. 

Measurement Result 
AcpiKsn+) w.r.t. D+ K-n+Tv+ (-1.6 ± 1.5 ± 0.9)% 
Acp{KsK+) w.r.t. D+ K--n+-n+ (±6.9 ± 6.0 ± 1.5)% 
AcpIksK+) w.r.t. D+ -» KsTi+ (±7.1 ± 6.1 ± 1.2)% 



our Monte Carlo. The Z)+ Kstt^ topology and recon- 
struction algorithm is substantially different from that of 
the K~TT^TT'^ and the two modes differ in how 

well the Monte Carlo matches to the data. For example 
there is a slight difference in how well the generated and 
accepted momentum distributions agree in each case. We 
use a technique modeled after the S-factor method used 
by the Particle Data Group ||l^ to evaluate the system- 
atic error. A scaled variance is calculated using the eight 
independent subsamples. The split sample systematic is 
defined as the difference between the scaled variance and 
the statistical variance when the former exceeds the lat- 
ter. Due to the smaller statistics in the _D+ KsK^ 
decay mode we can only form four independent subsam- 
ples. These are based on the run period in which the 
data were collected and on the normalized Ks mass. 

We evaluate systematic uncertainties due to the fitting 
procedure by calculating our results for various fit con- 
ditions, such as rebinning the histograms, changing the 
background shapes and in the case of KgK~^ also 

fitting the Dg peak. Since these different results are all a 
priori likely we use the resulting sample variance as a sys- 
tematic. The total systematic is calculated by adding the 
fit-variant systematic and the split-sample systematic in 
quadrature. For the Kstt^ measurements the sys- 

tematic has contributions from both the split-sample and 
fit-variant analyses. For the T{D+ K°tt+) / V{D+ 
K 



'7r± '7r± 



) measurement the contribution from the split- 
sample is 0.301% and from the fit-variant 0.098%. For the 
Acp{Ks'K^) measurement the split-sample contribution 
is 0.92% and that of the fit- variant is 0.13%. We find no 
systematic contribution to the KgK^ measure- 



ments from the split-sample technique, and therefore the 
fit- variant contributions are identical to the total system- 
atic error and are as shown in Tables || and |l[ Due to the 
lower statistics we did not split the KsK^ sample 

by momentum. Instead we treat the weighted average of 
two samples split by momentum as a fit variant. 

To conclude, we have searched for evidence of direct 
CP violation in the decays Kstt'^ and 

KsK~^ and measured their branching ratios relative to 
each other and to D'^ K^tt^tt^. Our relative branch- 
ing ratios are a considerable improvement over previous 
measurements. The CP asymmetries have not been pre- 
viously measured for these modes and are consistent with 
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